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The incomplete plasma dispersion function is a generalization of the plasma dispersion function in which 
the defining integral spans a semi-infinite, rather than infinite, domain. It is useful for describing the linear 
dielectric response and wave dispersion in non-Maxwellian plasmas when the distribution functions can be 
approximated as Maxwellian over finite, or semi-infinite, intervals in velocity phase-space. A ubiquitous 
example is the depleted Maxwellian electron distribution found near boundary sheaths or double layers, 
where the passing interval can be modeled as Maxwellian with a lower temperature than the trapped interval. 
The depleted Maxwellian is used as an example to demonstrate the utility of using the incomplete plasma 
dispersion function for calculating modifications to wave dispersion relations. 

PACS numbers: 52.35.Fp,52.25.Mq,52.40.Kh 



I. INTRODUCTION 

The incomplete plasma dispersion function, defined as 

Z {y, W ) = — / dt- , (1) 

V7T J v t-W 

for ^{w} > and as its analytic continuation for ^s{w} < 
0, was introduced by Franklin. 1 It is a generalization of 
the plasma dispersion functionpEI which corresponds to 
the limit v — > — oo: Z(w) = Z{— oc,w). Like the plasma 
dispersion function, it arises in the linear plasma di- 
electric function and corresponding wave dispersion rela- 
tions. However, the plasma dispersion function describes 
Maxwellian distributions, whereas Eq. describes non- 
Maxwellian distributions so long as they can be approx- 
imated as Maxwellian in finite, or semi-infinite, intervals 
of velocity phase-space. Each interval may have different 
characteristic densities, flow speeds and temperatures as- 
sociated with them. 

Situations where a plasma is far from equilibrium, but 
model distribution functions can be formulated in terms 
of piecewise Maxwellians are prevalent EHM1 An early ap- 
plication was in the study of waves near the boundaries 
of Q-machine plasmas. 8 12 There ions were modeled with 
a distribution consisting of only the high-energy tail of 
a Maxwellian. Electrons trapped by the confining ion 
sheaths were considered Maxwellian, whereas the pass- 
ing interval was depleted in density. These distributions 
are depicted by the dashed and sold lines, respectively, 
in Fig. [IJl 

This type of depleted electron distribution is com- 
monly found in the presence of potential barriers, such as 
sheaths near material walls or pr ob es as we n 
as double layers which provide an electrostatic barrier 
between plasmas with differing properties (density, tem- 
perature, etc. )J2SH!U Potential barriers create a trapped- 
passing boundary for one or more species. At the edge 
of an absorbing wall sheath, for example, the passing in- 
terval of the electron distribution will be empty. This 



depleted region can significantly alter the plasma dielec- 
tric response. For example, it has been shown to reduce 
the threshold for temperature-anisotropy-driven whistler 
instabilities near the boundary of magnetized plasmas 
Scattering can fill in the depleted region such that the 
density of the passing interval increases away from the 
potential barrier. Also, in some applications, secondary 
electron emissi on can generate an additional population 
of tail electrons EHEH see Fig.[l]3, which may also be mod- 
eled using Eq. 0. 

Other examples include electron sheaths, double 
sheaths and double layers that arise when a boundary 
is biased more positive than the plasma potentialP^LU 
In these situations the electron distribution can often 
be modeled as a depleted half-Maxwellian with a flow 
shift; see Fig. [IJl Phelps and Allen studied waves in 
the presence of a double sheath with electron emitting 
boundaries!^ They modeled the electron distribution as 
a background Maxwellian with an additional flowing half- 
Maxwellian component generated by wall emission. By 
modeling the electrons in this manner, they were able to 
accurately predict measured wave dispersion relations for 
high frequency waves. They noted that ions in these con- 
figurations can resemble Maxwellian distributions with a 
hole at low energies, as shown in Fig. 

The incomplete plasma dispersion function may also 
be useful for studying waves after non-linear process have 
modified a velocity distribution function, such as flatten- 
ing of a region of velocity-space due to a driven wave or 
instability,^ as shown in Fig. [l^. Another common ex- 
ample is the plateau distribution following the quasilinear 
evolution of a bump-on-tail distribution, as depicted in 
Fig. [T]i. These have been mea sure d for ions downstream 
of current-free double layersp32H Piecewise Maxwellian 
distributions have also been used to model ele ctrons up- 
stream of current-free double layers ) 22 ! 25 ! 26 * 28 ! as well as 
double layers that form in electronegative plasmas.^ 

The incomplete plasma dispersion function may also 
find use in magnetic reconnection research. A recent the- 
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passing boundary affects the bulk distribution (i.e., it is 
near the thermal speed), and the temperature charac- 
terizing the passing population is low enough. For in- 
stance, one significant modification to Langmuir waves 
is the absence of Landau damping for modes with phase 
velocity beyond the trapped-passing boundary because 
of the lack of resonant electrons. Similar undamped (or 
weakly damped) waves have previously been studied in 
the high frequency limit near a floating probe using a to- 
tal depletion model. The modeled wave properties were 
also confirmed using a probe and spectrum analyzerP 
Similarly, the dispersion relations developed in Sees. |III| 
and |IV| should be accessible to experimental measure- 
ment with the classical probe techniques. In addition, 
modern laser-induced fluorescence techniques may be ca- 
pable of measuring both the dispersion properties of these 
waves, as well as how they may influence the distribution 
functions in velocity-space pSMSl 



II. PROPERTIES 



FIG. 1. Example distributions functions for which the incom- 
plete plasma dispersion function will arise in the linear di- 
electric: (a) depleted Maxwellians with a flow-shift (dashed) 
and without (solid), (b) depleted Maxwellian with an addi- 
tional tail population from secondary electron emission, (c) 
Maxwellian distribution with a hole at low energy, (d) plateau 
distributions, (e) Maxwellian with a narrow flattened interval, 
and (f) a flat-top distribution. 



ory proposes that large scale parallel electric fields can 
both accelerate electrons and create a trapped-passing 
barrier that affects the electron distribution pSMU j n ^his 
model, the parallel electron distribution is Maxwellian 
in the tails and constant in the intermediate region, as 
shown in Fig. [1J. Notably, Eq. ([!]) will arise in the dis- 
persion relation of collisionless tearing modes, which can 
be a trigger for fast reconnection in magnetospheric plas- 
mas. 

The purpose of this paper is to provide several prop- 
erties of the incomplete plasma dispersion function that 
are useful for applying it to the linear dielectric response 
and wave dispersion relations in piecewise Maxwellian 
plasmas. These properties are discussed in Sec. [TTJ In 
Sees. |III| and |IV[ the utility of this function is demon- 
strated by applying it to the application of electrostatic 
waves near absorbing boundaries of unmagnetized plas- 
mas. In Sec. |III[ the depleted region of the electron dis- 
tribution is taken to be completely devoid of particles, 
whereas in Sec. |IV| the depleted region is modeled as a 
colder Maxwellian (returning the results of Sec. Ill in the 
limit of zero tail temperature). Waves in plasmas with 
gently deplet ed t ails ft^ as well as Lorentzian (or kappa) 
distribution J 41 * 42 1 have also been studied previously. 

Ion-acoustic and Langmuir wave dispersion relations 
are shown to be significantly modified when the trapped- 



A. Differential representation 

The incomplete plasma dispersion function, Eq. ([!]), 
can also be represented in the differential form 

_ 2 

dZ ^ 1 e~ v n / \ 

~ + 2wZ = — erfcM. (2) 

aw v — w 

Equation ([2| can be used to write higher order derivatives 
in terms of lower-order derivatives 

Z (n) = (rc-l)!e-" _ 2 r _ 1)z (n-2) wZ (n-l)l (3) 

for n > 2. Here Z^ = d n Z/dw n denotes the n th deriva- 
tive of Z with respect to w. 



B. Asymptotic expansions 

Asymptotic expansions can be obtained by applying 
the Plemelj formula 

2 i r°° e" £2 

Z(v,w)~iy/TrH(w-v)e- w \ dt , (4) 

Vtt J v t-w 

in which H is the Heaviside step function, then using 
a large or small argument Taylor expansion of the inte- 
grand. For \w\ <C 1, this yields 



Z(y, w) = iV^H(w - v)e~ w2 + J2 a ^ ( 5 ) 



in which 



r(-n/2, v 2 ) + 2a(n + 1, i/) 7 (-n/2, v 2 ) ^ 

a n = Q -^= • (6) 



3 



Here 

a(n,i/) = [l + (-lH(l-i//M)/4 (7) 

is a parameter that is unity if v < and n is even, but 
is zero otherwise, r(n, x) is the upper incomplete gamma 
function, and 7(71, x) is the lower incomplete gamma 
function™ The second term in Eq. (p| is taken to vanish 
for even n (due to a = 0), although t(— n/2, v 2 ) diverges 
for even n. The first four terms of the coefficient a n are 

a = E 1 {v 2 )/{2^), (8a) 
ai = e'^/ivVn) - erfc(i/), (8b) 
a 2 = e- iy2 /(2^ 2 ) - E l {v 2 )/{2^), (8c) 



l-2i/ 2 



«3 



^erfc(^), 



(8d) 



3^v 3 

where erfc is the complimentary error function and 
Ei(v 2 ) = T(0,^ 2 ) is the exponential integral™ 

Applying the same procedure in the asymptotic limit 
\w\ ^> 1 yields 



Z(y, w) ~ i(j\^7rH(w — v)e 
in which 



n=l 



b n /w n (9) 



r(n/2, ^ 2 ) + 2a(n + 1, v)j(n/2, v 2 ) 



and where 



0, > 0, 

a={l, 5S{w} = 0, 
2, 9{w;} < 0. 



(10) 



(11) 



Note that b n = a_„. The first five terms of the coefficient 
b„ are 



h = erfc(i/)/2, 

6 2 = e- 1 ' 2 /(2v / ^), 

6 3 = ^e- 1/2 /(2v / 7r) +erfc(i/)/4, 
b 4 = (l + v 2 )e-» 2 /(2^), 



3p/2 



-erfc(v). 



(12a) 
(12b) 
(12c) 
(12d) 

(12e) 



2^ 

C. Generalized functions 

Generalized functions of the form 

1 f°° t n e~ t2 

Z n {v,w) = —= dt , n>0, >0 (13) 

A/7T J u t-w 

can be related to derivatives of Z(y, w) via the generating 
function 



1 n 



1=0 

1 (_1 
7^ ^ 



(14) 



^-m-l)!(-l) m ff m (i,)e-- 



ra=0 



(^ — w) l ~ 



TABLE I. Values for the coefficients di(n) and c n (l). The top 
row gives values of n, and the left column values of / [e.g., 
<fe(5) = 20]. The Hermite polynomials are generated from 
Hi = H l n = Cn(l)x n [e.g., c 2 (4) = -48]. For Z = n the left 
value corresponds to c n (l) and the right value to di(n). 
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where H m denotes the m th Hermite polynomial, 
and di(n) denote the coefficients satisfying t n = 
2- n Y^= di(n)Hi(t)] see Tbl. |l| (from Tbl. 22.12 of 
Ref . I46j) . A derivation of Eq. jl4| is provided in appendix 
[A) The first few terms of Eq. (|14|) are 



Z (v,w) = Z, 
Z\{y,w) 



Z 2 (v,w) = 



Z' e~ v _ erfc(z/) 

Y + 2ypH{y-w) ~ 2 
2v(y-w) - 1 



2Z + Z" 
w[erfc(z/)/2 



V / 7r(^ — ^) 2 

■wZ]+e~ u2 /(2>/5F). 



(15a) 
^4l5b) 

(15c) 



D. Numerical evaluation 

The incomplete plasma dispersion function can be 
evaluated using similar techniques used to evaluate the 
plasma dispersion function. When the pole at t = w 
falls outside the integration interval (w < v), Eq. ([I]) can 
computed by direct integration. When the pole at t = w 
falls within the integrand (i.e., v < w < oo), a similar 
integral spanning — oo to v can be computed by direct 
integration, and the result subtracted from the plasma 
dispersion function to yield the incomplete plasma dis- 
persion function 



Z(v,w) = Z(w) 



i_ r 



dt- 



t 



(16) 



Efficient meth ods, such as continued fraction expansions 
(see Sec. |IIF), have been developed for evaluating Z(w). 



Figure [2] shows plots of the real and imaginary parts of 
the incomplete plasma dispersion function as a function 
of the complex argument w — x+iy for three values of the 
cutoff parameter v = —1,0 and 1. The complete plasma 
dispersion function is also shown as a dashed line for ref- 
erence. For cuts along the real line (y = 0) the real part 
of Z(y, w) has a pole at x = as shown in panel (a). As 
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FIG. 2. Real and imaginary components of Z(i/,w) in which 
w = x + iy. Solid lines correspond to various values of v\ 
v = — 1 (red), = (black) and v — 1 (blue). The dashed line 
shows the complete plasma dispersion function Z(— oc,w). 
(a) and (b) take y = and (c) and (d) take y = —0.5. 



w — >• ^, the function diverges logarithmically. This can 
be seen by considering the principal value term of Eq. Q , 
applying the substitution s = t — w and expanding the 
exponential term for s <C w in the dominant interval near 
s = 0. The leading term shows that the integral scales 
as — exp(— w 2 ) ln(|z/ — w\) as w — >• The imaginary part 
is equal to ^H(x — v) ex.-p(—x 2 ) for y = 0, as show in 
panel (b). For cuts along the line ?/ = — 0.5, the real part 
closely follows the complete plasma dispersion function 
for x > and no singularities arise, as shown in panel 
(c). The imaginary part is again discontinuous at x = z/, 
but has a more complicated behavior than the simple 
Gaussian shape found along the real line, as shown in 
panel (d). 

Figure [3] shows an altitude plot of the complete plasma 
dispersion function and incomplete plasma dispersion 
function for v = 1 in panels (a) and (b) respectively. 
The plot is obtained from the representation Z = \Z\e 10 
where blue lines are lines of constant \Z\ and red lines 
are lines of constant 0. In the incomplete case^ a branch 
cut arises at x = v creating a discontinuous boundary for 
constant \Z\ lines (see Refs. l47H49l for a discussion). The 
intersection of this and the real line locates a singular 
point in Z(y^ w) that is not present in the complete case. 



E. Two-pole approximation 

Because the plasma dispersion function is difficult to 
deal with analytically, it is often useful to have an ap- 
proximate algebraic expression. A method that has been 
shown to accurately capture both the real and imagi- 
nary parts of the conventional plasm a dispersion function 
is the two-pole approximation!^^ This method can be 
extended to the incomplete plasma dispersion function. 



The approach is to first write a Pade approximate of the 
form 



Z(y,w) 



1 + qiw + q 2 w 2 



(17) 



in which the coefficients are determined by matching 
the large and small argument asymptotic expansions of 
Eq. ([l7j) with those of Eq. ([!]). Using Eq. (9), the large 
argument limit gives pi/q2 = &i- Using Eq. (5), the small 
argument expansion of Z(y, w) is Z ~ Y^=o c nW n , where 



c = i\^7rH(x — u)-\- 



Ei{v 2 ) 
2v^ 



Cl 



-erfc(z/) (18) 



and 



C2 = —i\pKH(x — v) + 



2^ ' 



(19) 



Matching these with the small argument expansion of 
Eq. (17), the p coefficients are 



Po = c , Pi ~- 

and the q coefficients are 

c ci + c 2 b 1 



c c 2 ) 



qi 



-ci&i 



Q2 



Ci - c c 2 



Next, a two-pole form of Eq. (17), written as 
A 



Z(y,w) 



B 



(20) 



(21) 



(22) 



can be obtained from its partial fraction decomposition. 
This provides 



Qi 
2<z 2 



2q 2 



and 



Po +pic + 



B 



Po+PxC- 



(23) 



(24) 



Equation (22) is a two-pole approximation of the in- 
complete plasma dispersion function. In the complete 
limit {y —> — oo), the coefficients simplify considerably: 
A = 0.5- 1.289z, B = A* and c± = tO.514 - 1.032iPl 

A comparison of the two-pole approximation from 
Eq. p2|_ and numerical solutions of Eq. ([!]) are shown 
in Fig. |4^i. Here, the complete case (y = — oo) and an 
incomplete case (v = —0.1) are shown and the argument 
is taken to be real w = x. This Pade- type approxima- 
tion captures both the large and small argument limits, 
and smoothly connects them. However, the incomplete 
plasma dispersion function has a pole on the real line 
at w = v (see Fig. [2|i), which this approach does not 
capture. This leads to significant errors in the vicinity 



X 
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FIG. 3. Altitude plot of (a) Z(— oo,w) and (b) Z(l,w) where w = x + iy. 



of the pole. This technique has a limited range of va- 
lidity for the incomplete plasma dispersion function, and 
one must use it cautiously. On the other hand, as Fig. [4] 
shows, it can still be useful when the argument is far from 
the pole. A common situation where this approximation 
may be useful is when v < 0, but one is interested in 
waves for which wr > 0. 

Figure shows a measure of the relative error (Z — 
Z)/\Z\ between numerical solutions and two-pole approx- 
imations, again taking a real argument which here varies 
from to 20. The error grows as the cutoff parameter 
is increased from approximately 10% maximum in the 
complete case to approximately 30% for v = —0.1. 



F. Continued fraction expansion 

The two-pole approximation of the previous section 
was based on a Pade approximate written in a rational- 
fraction form 



N 

Z{y, w) ~ Pm(w) = ^2 A nW n I ^2 Bn%J ° n ( 25 ) The contmued fractions are of the form 

n=0 



M 



numerically evaluate the plasma dispersion function. In 
this section, a continued fraction expansion of Eq. 
and algorithm for evaluating it are calculated following 
the procedure presented in Ref. [52] (for an alternative 
method, see Ref.l53|)!^The result is also used to calculate 
a rational fraction representation. 

McCabe and Murphy^ have developed an algorithm 
for calculating a continued fraction expansion of a func- 
tion f(z) corresponding to power series expansions at two 
points. The expansions are assumed to have the form 



n=0 



and 



f(z) = -J2b n /w n 



(26) 



(27) 



with ao,ai 7^ 0. For our application, Eq. (26) will corre- 
spond to the small argument expansion and Eq. (27) to 
the large argument expansion of Eq. (fil). 

forrr 



n=0 



where N = 1 and M = 2. This approach can be extended 
to higher order to obtain a convergent approximation. 
Equivalently, the Pade representation can be written as 
a continued fraction. This is a common method used to 
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FIG. 4. (a) Comparison of the numeric ally evaluated Z and 
the two-pole approximation Z from Eq. ( 22 ) for the complete 
case (y — — oo) and an incomplete case (i/ = —0.1) taking 
a real argument w = x + Oi. The lines representing numeri- 
cal evaluation of the imaginary components lie on top of one 
another, (b) A measure of the relative error of the two-pole 
approximation [(Z — Z)/\Z\] for three values of the cutoff pa- 
rameter v — — oo, —0.4 and —0.1. Again, the argument is 
real and x : — »• 20. 



The coefficients n,- ^ and cL- ^ can be calculated from the 



recurrence relations^ 



(r) 
At) 



At) 



(r+1) 



Ar+1) 



At+1) 



n {r+1) d {r) 
n i+i a i 



(29a) 
(29b) 



for 



1,2,3, 



and 



values 



0, ±1, ±2, . . . using the initial 



-ao/^i, = — a r /a r _i for r > 1, 



(r) 



— b| r | /b| r - ii for r < — 1 and n{ = for all r. Solving 
Eqs. ( |29a| ) and (29b) provides an array of terms, called 
the n-d array. Once these are determined, the continued 
fraction of Eq. (28) can be evaluated to order N. 



McCabe and Murphy developed a convenient method 
to determine the convergents F\ , F2 , • • • , Fn , which also 



N=13 

° N=25 
Converged 




FIG. 5. Real (blue) and imaginary (red) components of the 
N th convergent of the continued fraction expansion of the 
incomplete plasma dispersion function Zn(v,w) for v — — 1 
and w = x — 0.5i. Solid lines show a converged solution 
obtained using the method of Sec. |II D| 



provides rational Pade approximates of the form of 
Eq. (25). The N th convergent (Fat), corresponding to 
the continued fraction of Eq. (28) keeping terms i = 



1,2,. 



, TV, can be written in the rational fraction form 
F N (w)=R N (w)/S N (w). (30) 



Here Rn and Sn can be determined from the recurrence 
relations^ 



Rn+i = (1 + d$ +1 w)R N + n K x } +1 wR N -i, (31a) 

S N+1 = (1 + dP +1 w)S N + nP +1 wS N -! (31b) 

for N > 1 and with the initial conditions Rq = 0, R\ = 
So = 1, and Si = 1 + d^w. 
Because the asymptotic expansions of Z(y,w) depend 
on H(wr — v\ it is convenient to use the representation 



Z(y,w) = i\/ttH{wr — v)e w + f{y,w) 



(32) 



and approximate f{v,w) with a continued fraction ex- 
pansion; /(z/, w) c± Fn(v, w). In this case, the coefficients 
of Eqs. (26) and ( [27] ) are simply the large and small ar- 
gument expansions from Eqs. ([6| and ( p~Q| ) (i.e., a n = a n 
and b n = 6 n ). Alternatively, the exponential term in 



Eq. (32) could be included in the expansion coefficients, 
but then the coefficients become functions of wr through 
the H(wr — v) term. 

Figure |5] shows convergents of Eq. (32) for N = 13 and 
25 obtained from the continued fraction method. Here 
v = — 1 and y = —0.5. Also shown is a converged numer- 
ical solution obtained using the method of Sec. |IID| The 
continued fraction solution converges quickly over much 
of the domain, with the exception of points in the vicin- 
ity of large gradients. Since it is very difficult to capture 



7 



regions with steep gradients using Pade approximates, 



we have found the method of Sec. II D to provide a more 
efficient means for numerical evaluation. The advantage 
arises from splitting the integral into a smooth part (the 
complete plasma dispersion function), which can be cal- 
culated efficiently with continued fractions, and a part 
with steep gradients, which can be evaluated by direct 
integration. 



III. WAVES IN THE PRESENCE OF A DEPLETED 
ELECTRON DISTRIBUTION 

Next, Eq. ([I]) is applied to calculate the dispersion 
relation of linear electrostatic waves near an absorbing 
boundary. The electron distribution in this scenario will 
have a trapped-passing boundary corresponding to the 
local electrostatic potential; see Fig. [6^1. Because the 
boundary is completely absorbing, the passing interval 
will nominally be empty. Collisions act to fill in the de- 
pleted interval, which populates as a function of distance 
from the boundary. In this section, we take the passing 



IV 



interval to be completely devoid of particles. In Sec 
we will generalize this to include a passing population 
with an independent temperature. 

The dispersion relation of linear electrostatic waves can 
be computed from the roots of the dielectric function 



1 



k 2 m s 



3 k • V v / S 

a v . 

uj — k • v 



(33) 



In this section, we take the electron distribution function 
to be 



fe = 



niH(v x - v c ) __ v 2 /v 2 



r 3/2 ? ,3 



(34) 



Tl 



see Fig. |6^l. In Eq. (34), n\ and T\ are parameters asso- 
ciated with the Maxwellian region of velocity space. The 
subscript 1 is used to denote this interval. In this case, 
the passing interval (interval 2) is empty so we need not 
consider it. We use this notation for consistency with the 
next section. The parameters rt\ and T\ can be related 
to the density, flow velocity, and temperature of the dis- 
tribution through the moment definitions: n s = J d 3 v f s , 
V S = J d 3 vvfJn s andT s = / d 3 v m s (v- V s ) 2 / s /(3n s ). 
Inserting Eq. (34) provides 



n e erfc(^i) 
th 2 ' 



and 



V T1 



= 1 



ri\ e 



x, 



2i/i V e 2 V 2 



3v^ 



3 vti 



(35) 



(36) 



(37) 




(c) 0.6 



FIG. 6. (a) Depleted Maxwellian distribution for z/i = 
v c /vti = -0.75 and four values of T 2 /Tn 0, 0.3, 0.6 and 
1. Also shown are the corresponding fluid moment parame- 
ters: (b) density, (c) flow speed, and (d) temper atur e. The 
case of total depletion (T2/T1 = 0) is used in Sec. 



Ill 



for the electron density, flow velocity and temperature, 
respectively. These are shown in Fig. [6] as a function of 
the cutoff parameter v\ — v c /vti- The coordinates are 
aligned so that v\ < 0. 

We consider only waves that propagate perpendicular 
to the boundary, such that k = kx. Taking the model 
distribution from Eq. (34) for electrons, and a flowing 



Maxwellian distribution for the ions, the dielectric func- 



tion from Eq. ( 33 ) reduces to 



k 2 v 2 Ti 



Z' 



k • Vi 



kv T i 



k 2 v\ x 



Z' 



U kvrij' 
(38) 

Next, numerical solutions and asymptotic approxima- 
tions will be used to study the dispersion relations of 
ion-acoustic and Langmuir frequency waves determined 
from roots of Eq. (38). We assume the flow is only in 
V 7 - = 



the parallel direction, = ViX, 
notation k • = kVi is used. 



and the abbreviated 



A. Ion-acoustic waves 

An approximate analytic dispersion relation for ion- 
acoustic waves is obtained from the large argument ex- 
pansion of the plasma dispersion function for the ion term 
(u — k • Vi)/(kvTi) ^> 1, and the small argument ex- 
pansion of the incomplete plasma dispersion function for 
the electron term uj/(kvTi) <C 1 using Eqs. ([2| and ([5|. 
Applying these asymptotic expansions to Eq. ( [38] ), the 
dielectric function reduces to 



e = l 



a u pi iSw 

k 2 \ 2 D1 ~ (uj - kVi) 2 + k 2 X 2 D1 



(39) 



LO r /kc s i 

FIG. 7. Cuts of the intersection of the real and imaginary 
parts of the electrostatic dielectric function with the [0, 0] 
plane. The real part is shown for k = 10 ~ 2 (red lines) and 
k = 10 2 (green lines). The imaginary part is the same for each 
k (blue lines) . The black lines show the dispersion relation of 
various modes as k varies from 10 -2 to 10 2 . 



in which w = oo/{kvTi), 
1 



a = 



erfc(^i) 



and 



5 y/nH(w R - vi) 



i (E x (y 2 ) e~ v2 



2 V V^ 1 



,2 / • 



(40) 



(41) 



Assuming that Sw <C l + /c 2 A| )1 , and that uor ^> 7, the 
dispersion relation obtained from the roots of Eq. ( 39 ) is 



IT iZi 



7rm e H(w R -v{) 
m V SMi {a + k 2 X 2 m f/ 2 



in which 



ur,± = kVi ± 



Zikc s i 



ni y/ a + k 2 X 2 m 



(42) 



(43) 



is the real frequency of the waves. The notation in 
Eqs. (42) and (43) assumes either the top or bottom set 



of ± symbols is used (not any permutation). The con- 
ventional ion-acoustic dispersion relation is returned in 
the limit n e /ni — » 1, a — >• 1, and H(wr — v) —> 1. 

Numerical dispersion relations have been obtained 
directly from Eq. (38). For ion- acoustic frequency 
waves, it is convenient to write Eq. (38) in terms of 



the dimensionless frequency and growth rate uj r jkc s \ 
and j/kcsi. Six free parameters then remain: 
k\Di, v\, Ti/Ti, Vi/c s i, m e /rrii and Z^. The dispersion 
relation for various roots at fixed k\r>\ are obtained from 
the points where both the real and imaginary parts of the 



10"' 



10 u 

k\ D] 




FIG. 8. The ion- acou stic dispersion relation computed numer- 
ically from Eq. (38) (solid lines) and from the approximate 
analytic formula from Eqs. (42) and (43) (dashed lines) for 
three values of the cutoff parameter v\ — —00, —1 and —0.1. 
The parameters T\/Ti — 100, Zi — 1, and mi/m e — 1836 
were used in each case. Panels (a) and (b) show the real and 
imaginary components for no flow V* = 0, while (c) and (d) 
correspond to sonic flow Vi = c s \. 



dielectric function, Eq. (38), vanish. Figure [7] shows cuts 
of each component in the[0, 0] plane. Here the red and 
blue lines denote the real and imaginary components for 
k\Di — 10 -2 . The other free parameters chosen in this 
plot are v\ = —00, T\/Ti = 10, Vijc s \ =0, rri z /rri e = 
1836 and Zi = 1. The upper set of circles denote the roots 
for this wavenumber, which each correspond to a differ- 
ent wave. Dispersion relations are obtained by stepping 
kXni and determining the root at each step. This proce- 
dure maps the black lines in Fig. [7] as kXjji varies from 
10 -2 to 10 2 . The green lines show the cut of the real com- 
ponent of i for kXpi = 10 2 and the imaginary component 
is essentially unchanged (blue lines). Figure [7] shows that 
these ion-acoustic dispersion relations lie along the zero- 
plane cuts of the imaginary component of e. 

Figure [8] shows example dispersion relations obtained 
with this method for three values of the cutoff parameter: 
v\ = v c /vti = — 00, —1 and —0.1. Here the frequencies 
are normalized to kc s \. Panels (a) and (b) take Vi = 
and look at the least damped root with uj r positive, while 
(c) and (d) take Vi = c s \ and look at the root where 
u r /kc s i is small, which is unstable for a range of kXm 
(there is also a higher frequency mode with uo r ~ 2kc s \ in 
this case, which is damped). The other parameters used 
are T\/Ti = 100, mi/m e = 1836 and Zi = 1. Solid lines 
represent the numerical solutions and dashed lines the 
analytic approximations from Eqs. (42) and (43). The 
curves for v\ = — 00 correspond to the conventional ion- 
acoustic wave in a Maxwellian plasma. The figure shows 
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that this is not significantly modified for v\ as large as 
— 1. Thus, modifications to the ion- acoustic dispersion 
relation due to the depleted electron interval near floating 
surfaces is negligible. For floating surfaces, v\ is typically 
a few times —1. Significant (~ 50%) corrections arise in 
the v\ = —0.1 case. Such a case is relevant near isolated 
boundaries, such as probes, biased so that the potential 
drop from the plasma to probe is less than the float- 
ing potential. Boundaries biased more positive than the 
plasma often result in electron distributions with v\ = 0. 

Figure |8| als o shows that the analytic approximations 
from Eqs. (|42[) and (43) are typically within the expected 



0(y/Ti/Ti) ~ 10% accuracy for k\]ji < 1. One excep- 
tion is the growth rate for the v\ = —0.1 case in panel 
(d), which is off by ~ 50% for small fcA^i. One inter- 
esting feature of the depleted interval is that it can have 
a greater affect on the growth rate than it does the real 
frequency; as shown in panels (c) and (d). For v\ = —1, 
depletion has a negligible affect on the real frequency, but 
destabilizes low-fc waves. However, the depletion pro- 
duces competing effects, as demonstrated by the result 
that increasing v\ further to —0.1 enhances the stabi- 
lizing effect of finite fcA^i. This may be because the 



electron-ion temperature ratio, see Eq. (37), is substan- 
tially reduced in this case. 
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FIG. 9. Contours of R{£} = (red) and %{i} = (blue) 
in the complex frequency plane (uj = ujr + £7) for Langmuir- 
frequency waves from Eq. (45) taking /cAdi = 0.2. 



B. Langmuir waves 



To obtain an approximate analytic dispersion relation 



for Langmuir waves (uj 



o), the ion term in Eq. ([38 



is 0(m e /rrii) and will be neglected. Taking the large 
argument expansion of the electron term from Eq. ([9|, 
Z' ~ [erfc(^)/2 + 5(w)]/w 2 , where S keeps terms up to 
fourth order terms, and w uses the lowest order solution 
of uj (w = oJpe/kvTi)- Using these in Eq. (38) provides 
the dispersion relation 



2 i 7 2 2 



2kV e 



J pe 



(44) 



in which 



pe 



eric(i/)uj p l 1 /2 is the plasma frequency 



based on the density moment and V e is the electron flow 
moment speed from Eq. (|36l). 



The ion term in Eq. (38) is also neglected for the 



numerical analysis of Langmuir waves, which are deter- 
mined from the roots of 



£ = 1 



UJ, 



Pi 



k 2 v\ x 



Z' 



"I, 



x 



kVTl J 



(45) 



Figu re [9| shows some interesting features of the zeros of 
Eq. (|45[) as the cutoff parameter increases from —00. In 
these plots, kXni = 0.2 was chosen. In the complete 
case (vi = —00) the two usual Langmuir roots are found 
uj = dzcopi. The dispersion of each of these is identical, ex- 
cept for the sign of the phase speed. For v\ = — 2 (at this 
kXoi), these two weakly-damped waves persist nearly un- 
affected by the depleted region in velocity-space. How- 
ever, several damped modes are strongly affected by the 



depletion, leading to an asymmetry in ±ujr that is not 
present in the complete case. As v\ is increased further 
to —1, this asymmetry extends to the usual Langmuir 
roots as well. Here \uJr-\ < but only slightly so 

for v\ = — 1. A more substantial asymmetry in the real 
frequency of the two oppositely directed modes is found 
as the cutoff parameter increases further. 



Figure 10 shows the Langmuir wave dispersion rela- 
tion for each of these modes at two values of the cutoff 
parameter: v\ = —00 and —0.1. The complete case is 
symmetric, so ujr+ = —ujr-, as is typical. As the cutoff 
parameter increases two effects set in: (1) there is a shift 
in the frequency at asymptotically low wavenumber, and 
(2) the two oppositely directed modes no longer have the 
same absolute frequency. Effect (1) is simply due to the 
ujpi normalization. Equation (44) shows that uj ~ ±uJ pe 
at asymptotically small wavenumbers, where uj pe is the 
electron plasma frequency based on the total electron 



density. For v\ 



-0.1, 



Juj p i \prQnx « l/y/2, 



which accounts for this shift. Effect (2) sets in at higher 
wavenumbers, which the asymptotic theory does not cap- 
ture. Here, the magnitude of the real frequency can be 
significantly smaller for the mode with phase velocity in 
the depleted region, in comparison to the oppositely di- 
rected mode. Because the electron density vanishes in 
this region of velocity phase- space, their are no resonant 
electrons, and hence no Landau damping. This mode is 
undamped, even for large fcA^i, in contrast to the strong 
damping that affects the mode with positive phase veloc- 
ity. This can be seen in Fig. [lOb. 
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Although the continuous case of Eq. (47) is considered 
here, the results can be generalized in a straight-forward 
way to electron-emitting boundaries by applying a model 
for n 2 that accounts for this additional population of elec- 
trons in the passing interval. 



Applying Eqs. (46) and (47) to the moment definitions 



of the fluid variables provides 



n e erfc(^i) 
n[ 2 



n 2 erfc(-z/ 2 ) 



Tlx 



and 

Te 



V T 1 



111 

2Ti 



m (i-T 2 /T 1 ) c _^ i 



(48) 



(49) 



VTl 



Vti 



-erfc(z^i)e 1 ' 1 



(50) 

for the electron density, flow velocity and temperature, 
respectively. These are shown in Fig. [6] as a fun ction o f 
the cutoff parameter v\ = v c /vti (note v 2 — v\ \/Ti/T 2 ) 
for a few values of the temperature ratio T2/T1. Putting 
Eq. (46) into (33), the linear plasma dielectric takes the 
form 



k2v Ti 



Z' 



k • Vi 



kv Ti 



(51) 
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FIG. 10. Langmuir wave dispersion relation calculated nu- 
merically from Eq. (45) (solid lines) and analytically from 
Eq. (44) (dashed lines)~(a) shows the real frequency for pos- 
itive phase velocity (gj+) and negative phase velocity (cj_) 
modes, (b) shows the growth rate for each mode. The damp- 
ing rate for uj±(isi — —00) and cj+(zyi = —0.1) lie on top of 
one another. 



Next, the ion-acoustic and Langmuir waves from the pre- 
vious section are revisited using Eq. (51), which includes 
a non-vanishing tail population. 



Ion-acoustic waves 



IV. WAVES IN THE PRESENCE OF A PARTIALLY 
DEPLETED ELECTRON DISTRIBUTION 

Scattering and ionization cause the passing interval of 
the electron distribution to be refilled as a function of 
distance from the boundary. Here, the discussion of the 
last section is generalized to include a passing population 
(labeled population 2), characterized by T 2 < T\. The 
electron distribution then takes the form 



fe 



i/v* 



'; 2 
j t1 



f n 2 exp(-vl/v^ 2 )/v T2 , 
\ niexp(-^/^ i )^ T1) 



< v c 



(46) 



where v c < 0. For absorbing boundaries f e is expected 
to be continuous, in which case 



To obtain an approximate analytic dispersion rela- 
tion, we apply to Eq. ( [51] ) the usual ordering: (uj — k • 
Vi)/kvTi ^> 1, oj/kvTi <C 1 and uj/kvT2 <^ 1- How- 
ever, we note at the outset that the approximate dis- 
persion relation obtained in this manner cannot be ex- 
pected to reduce to Eq. (42) in the T 2 — >> limit because 
this would violate the uJJkvT2 <C 1 assumption. Esti- 
mating uj/kvT2 ~ Csi/vt2, the uo/kvT2 ^> 1 ordering 
requires T 2 ^> T\m e jmi. Following the same procedure 
as Sec. |III A] provides the approximate ion-acoustic wave 
dispersion relation 



in which 



n e 7rm e 



e(w fl ) 



niV 8Mi (a + k 2 X 2 D1 ) 3 / 2 



(52) 
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FIG. 11. The ion- acoustic dispersion relation computed nu- 
merically from Eq. (51 ) (solid lines) and from the approximate 
analytic formula (dashed lines) for values of the temperature 
ratio: T2/T1 = 1, 0.6, 0.3 and 0. The parameters v\ — —0.1, 
Ti/Ti = 100, Zi = 1, and mi/m e — 1836 were used in each 
case. Panels (a) and (b) show the real and imaginary com- 
ponents for no flow V% — 0, while (c) and (d) correspond to 
sonic flow Vi = c s i. 
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FIG. 12. Contours of ^{e} = (red) and %{i} = (blue) 
in the complex frequency plane (uj = ujr + 27) for Langmuir- 
frequency waves from Eq. ( |57[) taking /cAdi = 0.2 and v\ — 
-0.5. Panel (d) solves Eq. fej for the T 2 = case. 



B. Langmuir waves 



is the real frequency. Here 

erfc(^i) T\ri2 erfc(— ^2) 



2 T 2 ni 2 



(54) 



and 



/ rp \ 3/2 

e = H(w R1 -v 1 )+H(v 2 -w R2 ) -M (55) 

The conventional ion-acoustic dispersion relation can be 
returned in the T 2 /Ti — >> 1 limit (in which case n e /ni, a 
and 9 1). 

Numerical solutions of the dispersion relation were ob- 
tained from Eq. (51) in terms of uj/kc s \ using the same 
method described in Sec. Ill A| Here seven free parame- 
ters must be set: kXoi, v\, T2/T1, Ti/Ti, Vi/c 3 i, m e /rrii 
and Zi. Figure |TT| shows the real frequency and growth 
rate as a function of kXpi for a few trapped-passing tem- 
perature ratios for stationary (Vi = 0) and sonically flow- 
ing (Vi = c s i) ion distributions. The other free parame- 
ters chosen were: v\ = —0.1, T\/Ti = 100, rrii/m e = 
1836 and Zi = 1. Figure [TT] shows a transition be- 
tween the conventional ion-acoustic wave (T 2 /Ti = 1) 
and the fully depleted wave (T 2 /Ti = 0) as T 2 /Ti de- 
creases. It also shows that the approximate formula from 
Eq. (52 ) typically captures the small kXpi behavior of the 



real freque ncy and growth rate to within the expected 
Q(yjT e jTi) ~ 10% accuracy. Equation (42) was used for 
the T 2 /Ti =0 case, since Eq. (52) is not expected to hold 
in this limit. 



An approximate analytic expression for Langmuir 
waves can be obtained from Eq. (51 ) using the same pro- 
cedure as Sec. |IIIB} Neglecting the ion term, and obtain- 
ing up to fourth order in the large argument expansion of 
the electron terms, this procedure provides the dispersion 
relation 



'■u pe + -k v T1 



(56) 



2kV e 



-kv T i v\ 



/TWc(-z/ 2 ) 



in which u pe is the electron plasma frequency based on 



the total electron density n e from Eq. (48) and V e is the 



flow- moment speed from Eq. (49). Equation (56) reduces 



to the conventional Langmuir wave dispersion relation in 
the limit T 2 /Ti -> 1 (V e -> 0) and to Eq. (|44| in the limit 
T 2 -> 0. 

As in Sec. |IIIB[ the ion term is dropped for the nu- 
merical analysis for the high-frequency Langmuir waves, 
so Eq. (51) reduces to 



k 2 v\ x 



T&2 

T 2 n 1 



[Z'(w 2 )-Z'(v 2 ,w 2 )] 



Figure 
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(57) 

shows how the zeros of and ^{^}, and 

the corresponding dispersion relations, change as the 
trapped-passing temperature ratio is varied. The pa- 
rameters v\ = —0.5 and k\r>\ = 0.2 have been cho- 
sen for this figure. The usual symmetric roots are ob- 
tained for the complete case (T 2 /Ti = 1). As the 
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FIG. 13. Langmuir wave dispersion relation calculated nu- 
merically from Eq. (57) (solid lines) and analytically from 
Eq. (56) (dashed lines) for v\ — —0.1 and three values of the 
temperature ratio: T2/T1 = 1, 0.2 and 0. (a) shows the real 
frequency for negative phase velocity modes, and (b) shows 
the growth rate for each mode. 



temperature ratio drops, this symmetry is broken and 
the density of damped roots increases substantially for 
ur/uj p i < uikXoi ~ —0.1. As the temperature ratio 
drops, the asymmetry between the oppositely directed 
modes increases, and the Landau damping decreases. 

Figure [13] shows the Langmuir wave dispersion rela- 
tion for v\ = —0.1 and three values of the temperature 
ratio T2/T1 = 1, 0.2 and 0. The general behavior is 

The 



similar to the totally depleted case from Fig. [10 
T2/T1 =0.2 line shows a transition between the two lim- 
iting cases. Figure [T3]o shows that the undamped modes 
extend to higher wavenumbers for cold, but nonzero, tail 
temperature. However, Landau damping does not vanish 
completely for finite tail temperature. The approximate 
formula from Eq. (56) accurately predict the wave fre- 



V. SUMMARY 

Plasmas, especially in the presence of potential bar- 
riers, can be far from equilibrium. Nevertheless, mod- 
els can often be constructed to describe the distribu- 
tion functions using Maxwellians in finite velocity-space 
intervals, where each interval is characterized by dif- 
ferent effective densities, temperatures and flow speeds. 
Whenever such models are applied to calculate the linear 
plasma dielectric response, the incomplete plasma dis- 
persion function, denned by Eq. 0, will arise. Several 
properties of this function, such as asymptotic expansions 
and other approximations, that are useful for applying it 
to calculate wave dispersion relations were reviewed in 
Sec. [TT1 In Sees. [Till and [iVl this function was used to 



develop quantitative dispersion relations for ion-acoustic 
and Langmuir waves in plasmas near absorbing bound- 
aries. The passing interval of the electron distribution 
is depleted in density compared to the trapped interval 
in this circumstance. Substantial modifications to the 
dispersion relations were shown to arise if the passing in- 
terval is sufficiently cold, and the trapped-passing bound- 
ary is not too far onto the tail of the distribution. These 
conditions are particularly interesting near probes biased 
near (or above) the plasma potential. 



Appendix A: Derivation of generating function for Z n 



An outline of the derivation of Eq. (14), which relates 



Eq. (13) to derivatives of Z{y,w), is provided in this ap- 
pendix. First, it is useful to apply n derivatives to Eq. 
to give 



dt 



(t-w) 



n+l ' 



Applying 



e~ l _ 1 (d n /dt n )e~ 
(t - w) n+1 ~ n\ t-w 

-n-1 



(Al) 



(A2) 



d 
dt 



(n-m- 1)! (d m /dt m )e- t 



m=0 



n\ (t - w) n ~ 

and the relation 

(d n /dt n )e- t2 = (-l) n e- t2 H n (t), 



Eq. (Al) can be rearranged to show 



1 



dt- 



t — w 



(-l) n Z (n) 



(A3) 



(A4) 



(-l) n ^ (n-m- l)\(-l) m H m (u)e-^ 



m—l 

Using the identity 



(y — w) n ~ 



quency at small wave number. 



1 n 



(A5) 



1=0 



13 



in Eq. ( 13 ) shows 



1 v— ^ 1 f 

1=0 v Jv 



dt 



t — w 



(A6) 



derivation of Eq. (14). 



Finally, putting Eq. (A4) into Eq. (A6) completes the 
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